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CLAIMS: 

1. A method for reducing impulse noise in a signal processing system, 
comprising the steps of: 

estimating parameters of an alpha-stable distribution to model impulse 
5 noise that corrupts data signals input into a transmission medium of the signal 
processing system; 

sampling signals from the transmission medium; said sampling step 
storing the sampled signals in a memory; the sampled signals having a noise 
component and a data component; and 

10 computing with a prediction filter an estimate of the data components of 

the sampled signals using the estimated parameters of the alpha-stable 
distribution. 



2. The method according to claim 1, wherein said sampling step is 
15 performed from a digital subscriber line (DSL). 

3. The method according to claim 2, further comprising the steps of: 

recording samples of impulse noise signals {jci, jc 2 , ... , x n ] transmitted over 
the transmission medium of the signal processing system; the samples of 
impulse noise signals recorded by said recording step having no data 
20 component; and 
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using the recorded samples of impulse noise to perform said estimating 

step. 

4. The method according to claim 1, wherein said sampling step is 
performed using an imaging device. 

5 5. The method according to claim 4, further comprising the steps of: 

recording a first, a second, and a third image with the imaging device; 

computing a difference between the first and the second image to define a 
centro-symmetrized difference image; 

estimating a characteristic exponent of an alpha-stable distribution using 
10 the centro-symmetrized difference image; and 

computing a sample of impulse noise signals by applying a centralizing 
transformation to the first, the second, and the third images. 

6. The method according to claim 1 , further comprising the steps of: 

generating a synthetic noise signal with a random number generator that 
15 uses the parameters of the alpha-stable distribution as input parameters to the 
random number generator; and 

subtracting the synthetic noise signal from the sampled signal to produce 
a modified sampled signal; the modified sampled signal having a noise 
component with a symmetric distribution. 

20 7. The method according to claim 1 , further comprising the steps of: 
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estimating a characteristic exponent of the alpha-stable distribution that 
models impulse noise corrupting the sampled signal; and 

optimizing model coefficients of the prediction filter using the estimated 
characteristic exponent of the alpha-stable distribution. 

5 8. The method according to claim 7, wherein said optimizing step 

minimizes a p^-power error criterion to optimize the model coefficients of the 
prediction filter. 

9. The method according to claim 8, wherein said optimizing step defines 
the p^-power error criterion using the characteristic exponent of the alpha-stable 

10 distribution. 

10. The method according to claim 8, wherein said optimizing step 
minimizes the p^-power error criterion by performing the following computation 
iteratively: 

15 where c are model coefficients, w is a diagonal weight matrix, and x t is 

an observed signal block of l samples, where: 





a 














c= 


b 

c 


, such that a = 




, b = 


2 


, and c = 


6 



x, = 



x[txL] 
x[f xL + L-1] 



, where t = 0, 1, 2, 3, 
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X ea = [x {l) X w X 0) \, such that 



x (,) = 



x{txL] 
x[txL+l] 



0 

xfrxL] 



0 
0 



x[txL+L-i\ x[txL+L-2] ••• 4fxL+L-AT| 



X (2) = 



* 2 [rxL] 
x 2 [fxL-H] 



0 

4fxL + l]4fxL] 



x 2 [txL + N] x[txL + N]x[txL + N-l] 



X (3) = 



jc 3 [rxL] 
* 3 [/xL + l] 



0 

x 2 [txL + l]x[txL] 



x 3 [txL + N] x 2 [txL + N]x[txL + N-l] 



0 
0 

x 2 [txL] 



x 2 [txL + L-l] x[txL+L-l]x[txL+L-2] ••• x 2 [txL + L-N] 



0 
0 

x 3 [?xL] 



x 3 [*xL + L-l] x 2 [txL + L-l]x[txL + L-2] ••• x 3 [rxL + L-JV] 



1 1 . The method according to claim 8, wherein said step of estimating the 
characteristic exponent comprises the steps of: 

computing a moment of the alpha-stable distribution; and 

estimating the characteristic exponent using the computed moment. 

1 2. The method according to claim 1 1 , wherein said step of computing the 
moment of the alpha-stable distribution comprises the step of computing an 

absolute fractional lower order moment A p according to the following equation: 
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A p =— ]jT \x k \ p , where p is the order of the moment and X k is one sample 

from a sequence of noise samples {Xi, X 2 , ... , X n }. 

13. The method according to claim 12, wherein said step of estimating the 
characteristic exponent comprises the step of estimating the characteristic 

exponent by computing Z = logA* -logA^ using the computed absolute 
fractional lower order moment S p by performing the steps of: 

partitioning the sequence of noise samples into two parts U and V, with 
each part containing data samples Ui, U2, U3, ... and V1, V 2) V 3 , respectively; 

computing the moment a* by summing noise samples in the sequence of 
noise samples as: 

Xi=Ui+Vi, X 2 =U2+V 2 , X3=U 3 +V 3 , and 

computing the moment A£by concatenating noise samples in the 
sequence of noise samples as: 

Yi=Ui, YaFVi, Y 3 =U 2 , Y 4 =V 2 , Y 5 =U 3 , Y 6 =V 3 

14. The method according to claim 13, wherein said estimating step 

estimates the characteristic exponent to be equal to ***° g2 when 

Z 

Z < plog2 and Z > /?log2 , where p is the order of the moment. 

<*min 2 
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1 5. The method according to claim 1 1 , wherein said step of computing the 
moment of the alpha-stable distribution comprises the step of computing an 
absolute logarithmic moment using the following equation: 

1 n 

L, =— ^ \og\x k \, where X k is one sample from a sequence of noise 

5 samples {X^ X 2) ... , X n }. 

16. The method according to claim 15, wherein said estimating step 
estimates the characteristic exponent by computing Z = Lf - L[ using the 
computed absolute logarithmic moment , by performing the steps of: 

partitioning the sequence of noise samples into two parts U and V, with 
10 each part containing data samples Ui, U 2) U3, ... and V1, V 2 , V 3| respectively; 

computing the moment £f by summing noise samples in the sequence of 
noise samples as: 

Xi=Ui+V 1( X 2 =U 2 +V 2 , X3=U 3 +V 3 , and 

computing the moment z£ by concatenating noise samples in the 
15 sequence of noise samples as: 

Y 1= Ui, Y 2 =V 1( Y 3 =U 2 , Y 4 =V 2) Y 5 =U 3 , Y 6 =V 3 , .... 

17. The method according to claim 16, wherein said estimating step 
estimates the characteristic exponent to be equal to 



when Z < and Z > 

"min 2 
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18. The method according to claim 16, wherein said computing step 
computes a signed logarithmic moment using the following equation: 



19. The method according to claim 18, further comprising the step of 
estimating a symmetry parameter of an alpha-stable distribution by computing 
tan((A/L,)a^/2)) 



20. The method according to claim 1, wherein said computing step 
extrapolates with the prediction filter an estimate of the data components of the 
sampled signals using the estimated parameters of the alpha-stable distribution. 

21. The method according to claim 1, wherein said computing step 
interpolates with the prediction filter an estimate of the data components of the 
sampled signals using the estimated parameters of the alpha-stable distribution. 

22. The method according to claim 1, further comprising the step of 
computing a moment of the alpha-stable distribution. 

23. The method according to claim 22, wherein said step of computing a 
moment comprises the step of computing one of: 

* 1 n 

an absolute fractional lower order moment A p \x k \ p , 



A = »X si 9 n ^) lo g|^l- 




tan(a^/2) 



a signed fractional lower order moment S p =— Y s\gn(X k )\x k 
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J n 

a signed logarithmic moment A = -^sign(X Jt )log|x Jt | l 

« jt=i 



1 n 

a first absolute logarithmic moment ^ = — log|x t | , 



a 



second absolute logarithmic moment L 2 = loglx*]-/^) 2 , 



a third absolute logarithmic moment L3 = (]£ loglxJ-Li) 3 , 

n-l t=1 



1 n 

an empirical characteristic function moment @(p) = —^e JpXk , 



a first upper extreme value moment Y x = Y.Kk , 

n/r * =1 



1 C 

first lower extreme value moment Y x = Y,£* > 

— n/rf^ 



a second upper extreme value moment Y 2 = -^i) 2 » and 

(n/r)-lS 



a second lower extreme value moment Y 2 = Y (K k -y,) 2 , 

— (n/r)-l^ — 



where: 

X* is one sample from a sequence of noise samples {Xi, X 2 

X„}, 

p is the order of the moment, 

Kk = max {log log X^., )+2 ,-~ log X^.,^, } , 
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K k =max{log -X^plog -X r(k . 1)+2 ,-- log -X^,.,}, and 

ris a block length of data. 

24. The method according to claim 23, further comprising the step of 
computing a characteristic exponent a of the alpha-stable distribution. 

5 25. The method according to claim 24, wherein said step of computing the 

characteristic exponent a of the alpha-stable distribution comprises the step of 
computing one of: 

a first auxiliary variable Z = log A* - log A Y p , 
p-; a second auxiliary variable Z = Lf - Lf , and 

# / ~ \-l/3 

io a third auxiliary variable Z = 1 - , and 

W 1.2020569 I 

V J 



a fourth auxiliary variable Z = 



2V6 



26. The method according to claim 24, further comprising the step of 
computing a symmetry parameter p of the alpha-stable distribution. 

27. The method according to claim 26, wherein said step of computing the 
15 symmetry parameter p of the alpha-stable distribution comprises the step of 

computing one of: 

a first estimate of the symmetry parameter B = tan ^ A/ ^ ) °^ 2 )) and 
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a second estimate of the symmetry parameter p - 1 - 



exp(a(Yi-Y x )) 



28. The method according to claim 24, further comprising the step of 
computing a dispersion y of the alpha-stable distribution. 

29. The method according to claim 28, wherein said step of computing the 
5 dispersion y of the alpha-stable distribution comprises the step of computing 



( r(l-p)cos(px/2) Y ' P 
J^il-p/a) cos{p6l a) 



cos0 1 , where 6 = arctan(/J tan(cr;z / 2)) , 



30. The method according to claim 24, further comprising the step of 
computing a location parameter 8 of the alpha-stable distribution. 

31 . The method according to claim 30, wherein said step of computing the 
10 location parameter 8 of the alpha-stable distribution comprises the step of 

computing a fractile / of the sequence of noise, where 

1 6 

f = , 6 = arctan(/? tan(OTr / 2)) . 

2 7ta 

32. The method according to claim 1 , further comprising the steps of: 
computing a characteristic exponent a of the alpha-stable distribution 

15 computing a location parameter 5 of the alpha-stable distribution, 

computing a symmetry parameter (3 of the alpha-stable distribution, and 

computing a dispersion yo\ the alpha-stable distribution. 
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33. The method according to claim 32, wherein said step of computing the 
characteristic exponent a and said step of computing the dispersion y of the 

alpha-stable distribution comprises the step of computing the following matrix 
equation 



j=m,k=m j= m,k =m 

Xlog(f,.)log(r t )C yt XlogC^C., 

;'=1.*=1 ;=U=1 

j-m,k-m j-mjc—m 

2ioge t )c, t £c„ 

y=U=l >=U=1 



a 
log y 



j=m t k-m 

j=m,k=m 



where 



{*i»*2»"*»'m} ' s some set of positive real numbers, 

C equals W 1 where W is an mxm matrix with elements 

w jk ={exp(f[tf +»f -rfff])-i}/{y a *7'f> 

and a\?*are existing estimates of the characteristic exponent and 
dispersion, and 
where y k = log(- log(Re <p(t k ))) for sequence Y for k = 1, • • • , m , and 

Y is a centro-symmetrized version of the sequence of noise 



samples. 



34. The method according to claim 33, wherein said step of computing the 
location parameter 8 and said step of computing the symmetry parameter p of 

the alpha-stable distribution comprises the step of computing the following matrix 
equation: 
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j=m,k~m j=m,k=m 

y=l,Jt=i ;=U=1 

j-m,k~m j=m,k=m 

* l<rc jk ^c ik 

j=\,k=\ j=\,k=\ 



j=m,k=m 

y"=U=l 
j=m,k=m 



where k - /tanf— 1 , 
V 2 J 



C equals W" 1 where W is the mxm matrix with elements 



1 



2{R 2 s +I 2 s ){R 2 r +I l T )t s t r 



{R s+r (I s I r -RM + Rs-rVJr+RsRr) 



+ I s+r (I s R r +R S Ir) + Is-r(IsK -RJ r )-*R s I s R r I r } 



where 



R s =Re(<p(t s ))J s = Im(?(f s )),R s+r = Re(?(r, +t r )\ 



where <p(t) = exp(-y\t\ a (1 + jfi'tan*f) + jtS') for existing estimates /J\<5' 
of the skew and location parameters, and 

0) k = -lmlog($(t k ))/ t k for the sequence of noise samples. 

35. An apparatus for reducing impulse noise in a signal processing 
system, comprising: 

a parameter estimation module for estimating parameters of an alpha- 
stable distribution; the alpha-stable distribution modeling impulse noise that 
corrupts data signals input into a transmission medium of the signal processing 
system; 
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a memory for accumulating sampled signals output from the transmission 
medium; the sampled signals having a noise component and a data component; 
and 

a signal estimation module for computing an estimate of the data 
5 component of the sampled signals output from the transmission medium using 
the estimated parameters of the alpha-stable distribution. 



36. The apparatus according to claim 35, wherein said signal estimation 
module further comprises: 

10 a prediction filter for estimating, using model coefficients, the data 

components of the sampled signals output from the transmission medium; the 
estimated data components corresponding to an estimation of the data signals 
input into the transmission medium; and 

a coefficient optimization module for optimizing the model coefficients of 
15 the prediction filter using one of the estimated parameters of the alpha-stable 
distribution received from said parameter estimation module. 

37. The apparatus according to claim 36, wherein the dependence of the 
prediction filter output on the model coefficients is linear. 

38. The apparatus according to claim 37, wherein said prediction filter is a 
20 Volterra filter. 
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39. The apparatus according to claim 36, wherein the dependence of the 
prediction filter output on the data components of the sampled signals input into 
said prediction filter is non-linear. 

40. The apparatus according to claim 36, wherein said parameter 
5 estimation module further comprises: 

means for computing moments of the alpha-stable distribution; and 

means for estimating the parameters of the alpha-stable distribution with 
the computed moments. 

41. The apparatus according to claim 40, wherein said parameters 
10 estimation module further comprises means for transforming the alpha-stable 

distribution to obtain deskewed alpha-stable random variables. 

42. The apparatus according to claim 40, wherein said parameters 
estimation module further comprises means for transforming the alpha-stable 
distribution to obtain centralized alpha-stable random variables. 

15 43. The apparatus according to claim 36, wherein said parameter 

estimation module adaptively estimates the parameters of the alpha-stable 
distribution. 

44. The apparatus according to claim 36, wherein said coefficient 
optimization module uses the parameters of the alpha-stable distribution to 
20 specify a minimum dispersion error criterion for determining the model 
coefficients of said prediction filter. 
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45. The apparatus according to claim 44, wherein said coefficient 
optimization module minimizes a cost function defined by the minimum 
dispersion error criterion to optimize the model coefficients of said prediction 
filter. 

46. The apparatus according to claim 45, wherein said coefficient 
optimization module minimizes a p^-power error criterion to optimize the model 
coefficients of said prediction filter. 

47. The apparatus according to claim 36, wherein transmission medium 
into which the data signals are input is a twisted pair. 

48. The apparatus according to claim 35, wherein the signal processing 
system operates in a digital subscriber line (DSL). 

49. The apparatus according to claim 35, wherein the signal processing 
system operates in an imaging device. 
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